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We present dielectric spectroscopy data obtained for gallium-doped Cdo.ggMno.oiTeiGa mixed crystals which 
exhibit a very special case of the two-power-law relaxation pattern with the high-frequency power-law exponent 
equal to 1 . We explain this behavior, which cannot be fitted by none of the well-known empirical relaxation 
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I. INTRODUCTION 

No doubt, the linear response theory ^ plays a great role 
in physics, being a theoretical background and fundamental 
accomplishments of statistical physics. Its conventional for- 
mulation is based on two assumptions |j2tl: 1) The time evo- 
lution of the system variables is governed by Hamiltonian op- 
erators; 2) The external perturbation arising from an initial 
disturbance makes the system depart from canonical equilib- 
rium weakly, and the linear response function is expressed in 
terms of the derivative of a stationary correlation function (of 
dipole orientation in the theory of relaxation, for example). 
This gives a rigorous and efficient approach to the temporal 
description of some Hamiltonian systems Jst]. However, the 
foundations are crashed on many-body kernel in dynamics of 
complex systems. The main cause is that the macroscopic 
evolution of such systems cannot be attributed to any particu- 
lar object chosen from those forming the systems. Any form 
of many-body interactions in the complex systems should be 
introduced through cross-correlations between their different 
objects, yielding hence insurmountable mathematical prob- 
lems. Therefore, it is very difficult, if possible generally, to 
describe the time evolution of many-body systems by Hamil- 
tonian operators, and any stationary correlation function is not 



available at all. Since, in general, the dynamical processes in 
complex systems are strictly stochastic in nature, they should 
be analyzed in a corresponding manner 

For description of relaxation and transport properties in 
such systems as glasses, liquid crystals, polymers, etc., the 
continuous-time random walk (CTRW) processes are one of 
the most useful mathematical tools. Despite their long history, 
started with the brilliant Montroll- Weiss idea Ij, the CTRWs 
are still far from their full exploration. However, their connec- 
tion with anomalous diffusion has been already recognized 
(see e.g., isj]). Recently, a progress in understanding of this 
mathematical tool ia-UOll stimulated new developments in dif- 
fusive scenarios of the non-exponential relaxation phenomena 

EMI. 
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In this paper we present experimental data which con- 
firm that the non-exponential relaxation behavior, in fact, 
is governed by compound counting processes strictly con- 
nected with clustered CTRWs |@]. In Sec. [Ill we study 
dielectric spectroscopy data measured for gallium-doped 
Cdo.ggMno.oiTeiGa semiconducting mixed crystals possess- 
ing deep, metastable defects. We observe that this material 
exhibits such a relaxation pattern which cannot be fitted with 
any of the well-known empirical relaxation functions. Hence, 
we apply a novel relaxation law being a modification of the 
result derived recently in ifisll . In Sec. |lll] we propose an 
anomalous diffusion scenario, based on the notion of a com- 
pound counting process, by means of which the observed re- 
laxation behavior may be explained. We end with conclusions 
in SecHVl 



2 



II. EXPERIMENT 

Dielectric spectroscopy studies carried out on various phys- 
ical systems revealed that a wide class of materials follows 
the anomalous relaxation mechanism ifTill represented by low- 
and high-frequency fractional power-law dependences of the 
imaginary part of the complex dielectric permittivity 

e*{uj)=e'{uj)~ie"{uj): 



(1) 



where Up denotes the loss peak frequency and < 7ti, n < 1. 
Depending on a mutual relation between the power-law ex- 
ponents two different types of relaxation responses can be 
distinguished. The relaxation response is called typical when 
the power-law exponents satisfy relation m > 1 — 7i. In or- 
der to interpret this type of experimental data the well-know 
Havriliak-Negami (HN) function IHlil 



[1 + [iuj/ujpY 



< a, 7 < 1 (2) 



is used: e* (w) = (eo - Eoo) ^hn* (w) + £oo, where eo is 
the static permittivity and represents the asymptotic value 
of the dielectric permittivity at high frequencies. For the HN 
function ^ the two-power-law property ([T]l is fulfilled with 
the power-law exponents m = a and 1 — n = 70;. For a long 
time period the HN function with extended parameters' range 
< a, Q!7 < 1 was also used to fit the less typical relaxation 
data for which the power-law exponents yield the opposite in- 
equality m < 1 — n. Unfortunately, none of the known relax- 
ation models lfl2l[T6ll can justify the values 7 > 1 appearing 
in the extended range of the power-law exponents. Recent 
progress in stochastic modeling of relaxation processes has 
resulted in derivation of a new relaxation pattern lfT2l [Tsll un- 
derlying the less typical responses: 



ifTsIl the relaxation response of Cdo.g9Mno.oiTc:Ga is mainly 
influenced by presence of deep, metastable traps within the 
band gap of the investigated Au-Cdo.99Mno.01Te Schottky 
junction. For purpose of the present study we analyzed the 
frequency-domain response of two samples of the same x = 
0.01 manganese (Mn) content labeled as sample 1 and sample 
2, respectively. Both the samples posses the same net donor 
concentration of approximately 10^ '' cm^'"* estimated from the 
capacitance-voltage measurements. Description of the sam- 
ples preparation can be found in details elsewhere |fl9ll. Gold 
Schottky contacts were thermally evaporated on the front side 
of the samples. Measurements were performed at zero bias us- 
ing Novocontrol impedance analyzer. Applied ac probe signal 
amplitude was equal to 10 mV. 

In Figure [U the normalized imaginary part of the dielectric 
permittivity for two samples, investigated in a broad temper- 
ature range, is presented. It is clear from the plot that both 
the samples exhibit the less typical, two-power-law relaxation 
pattern m < 1 — n, however, with different values of the low- 
frequency power-law exponent m. The change in values of 
this exponent, in the samples of the same manganese content 
and gallium concentration, may be associated with locally dif- 
ferent surroundings of DX centers contributing to the effective 
relaxation response of the investigated sample. 
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0<a,7<l. (3) 



This function exhibits the two-power-law property ([T]i with 
the power-law exponents m = and 1 — n = a. It not 
only properly describes the less typical class of relaxation re- 
sponses but also relates the experimentally observed power- 
law properties with random characteristics of the investigated 
system. Let us notice that both formulas ^ and ([3j can be 
used as fitting functions also with parameters a = 1 and/or 
7 = 1. In such cases they result from slightly different diffu- 
sion scenarios and may exhibit one fractional power-law only 
or even none. 

The less typical relaxation pattern is observed in semi- 
conducting mixed crystals of Cdo.99Mno.oiTe:Ga possess- 
ing deep metastable defects - the so called DX centers. It 
is widely accepted that depending on the position of gallium 
(Ga) dopants in the CdMnTe lattice shallow donor states or 
deep metastable traps may be formed iHtIi . As shown in Ref . 
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FIG. 1 : Normalized imaginary part of dielectric permittivity obtained 
for two samples of Cdo.ggMno.oiTeiGa. Both the samples exhibit 
two-power-law relaxation behavior with power-law exponents satis- 
fying relaxation m < 1 — n. 
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FIG. 2: Sample real and imaginary part of tlie permittivity data as a 
function of frequency obtained for tfie two investigated samples of 
Cdo.gtjMno.oiTeiGa at various temperatures. Solid lines represent 
function 
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FIG. 3: The values of fitting parameters obtained for two samples of 
Cdo.99Mno.oiTe:Ga obtained at various temperatures. Sample I : 
a = 1.00 ±0.01 and7 = 0.63 ±0.01; Sample 2 : a = 1.00 ±0.01 
and 7 = 0.53 ±0.01. 



In Figure |2] the experimental data fitted by means of for- 
mula Q are presented. It can be observed that this function 
perfectly covers the experimental data points. Values of the 
fitting parameters are collected in Figure |3] Both a and 7 
values are temperature independent. It should be pointed out 
that in case of both analyzed samples a parameter remains 
the same, whereas value of 7 is significantly different depend- 
ing on the sample considered. Since a is approximately equal 
to 1, we observe lack of the high-frequency fractional power- 
law behavior Such a relaxation pattern, being a very special 
case of ([T]), suggests a non-standard diffusion scenario under- 
lying the experimental result. 



III. DIFFUSION SCENARIO 

The CTRW process R{t) determines the total distance 
reached by a random walker until time t. It is characterized by 
a sequence of independent and identically distributed (i.i.d.) 
spatio-temporal random steps {Ri,Ti),i > 1. If we assume 
stochastic independence between jumps Ri and waiting times 
Ti, we get a decoupled random walk; otherwise we deal with 
a coupled CTRW. The distance reached by the walker at time 
t is given by the following sum 



i=l 



(4) 



where = max{n : X^'iLi — 0' counts the performed 
steps. 

Theoretical studies of the relaxation phenomenon in the 
above framework are based on the idea of an excitation un- 
dergoing (anomalous, in general) diffusion in the system un- 
der consideration isj]. The relaxation function (j>{t) is then 
defined by the inverse Fourier or the Laplace transform of the 
diffusion front R{t) sa R[t / tq) / f (tq) , where the dimension- 
less rescaling parameter tq ~ and /(tq) is appropriately 
chosen renormalization function. The diffusion front R{t) 
approximates a position at time t of the walker performing 
rescaled spatio-temporal steps (Ri/ f{To),ToTi). The char- 
acteristics of the relaxation process are related to the prop- 
erties of the diffusion front resulting from assumptions im- 
posed on the spatio-temporal steps of the random walk. For 
example, the decoupled CTRW with power-law waiting-time 
distributions (i.e., with the random variables Ti satisfying 
Prob(T:i > t) (t/to)"" as i 00 with some < a < 1 
and to > 0) leads to the Cole-Cole relaxation ll20ll . How- 
ever, the frequency-domain Cole-Cole relaxation with the cor- 
responding time-domain Mittag-Leffler pattern is only one of 
the cases measured in various experiments with complex me- 
dia, and derivation of those more general patterns requires 
considering diffusion scenarios based on a compound cou- 
pled CTRW representation. It should be pointed out that the 
simple coupling of type Ri ^ Tf (with positive power ex- 
ponent p) does not lead behind the Cole-Cole relaxation ll2lll . 
In contrast, introducing a dependence between the jumps and 
waiting times by a random clustering procedure we can ob- 
tain another empirical relaxation laws Uke the Cole-Davidson 
or Havriliak-Negami patterns lfT2i [T3II22I1 . Below, we present 
the diffusion scenario which leads directly to the results dis- 
cussed in the preceding section. 

Let Alj be a sequence of i.i.d. positive integer-valued ran- 
dom variables independent of the pairs {Ri ,Ti). Next, assume 
that the jumps and waiting times are assembled into clusters 
of random sizes Mi, A/2, • ■ •• This assumption allows one to 
transform the sequence of spatio-temporal steps {Ri, Ti) into 
a new sequence {Rj , Tj) of random sums 
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Mi+...+Mj 

iR,,f,) = J2 (^-^«)' J ^2. (5) 

i=Mi + ... + AIj-i + l 

Then the position R^^ (t) of the walker is determined by 
{Rj , Tj ) and, in accordance with the general formula (|4|i, it 
is given by 

= (6) 

where D{t) — max{n : X]j=i — 0- The dependence 
between the jumps Rj and the waiting times Tj of the coupled 
CTRW process is determined by the distribution of the 

cluster sizes Mj 

In the simple case when the waiting times are represented 
by equal intervals in time, i.e., Ti = At , we have 

[t/At] 

R{t)= ^ i?,, (7) 

i=l 

where [-J denotes the integer part, while the clustering pro- 
cedure yields Tj = MjAt for j > 1, and the coupled 
process R'^^t) in (|6]l takes an equivalent form 

R^'{t)^ J2 R^- (8) 
1=1 

Here U'^^{v{i)) is a compound counting process obtained 
from U^'{n) = Y^^^ f,/At = Y.]Li Mj, and v{t) = 

max{n : I]"=i < = max{ri : U^\n) < t/At}. 
Observe that formula dHJ is an analog of (|7]i with the com- 
pound counting process U'^^ {v{t)) substituting the determin- 
istic number [t/A^J of performed jumps Ri. The counting 
process U^^ {i>{t)) is always less than [t/At], and it is hence 
a special case of the undershooting compound counting pro- 
cess lfT2ll . It is also a clear signature of the spatio-temporal 
coupling provided by the clustering procedure (|5]l. 

The idea of compound counting processes in CTRW ap- 
proach is not new in physics. The resulting CTRW processes 
were examined in the context of the rareness hypothesis in the 
fractal-time random walk models (see, e.g. I23ll24ll ). In gen- 
eral, the compound counting process cumulates random num- 
ber of random events. Physical situations where the relevance 
of this scheme holds are numerous. For instance, take into 
account the energy release of individual earthquakes in geo- 
physics, the random magnitude of claims' sequence in insur- 
ance risk theory or random water inputs flowing into a dam in 
hydrology where summing the individual contributions yields 
the total amount of the studied physical magnitude over cer- 
tain time intervals. 

The diffusion front R^^ (t) related to (|8]l takes the subordi- 
nated form 

R'\t) ^ X{Z{t)) , (9) 




2^4 6 2 \ A 6 




FIG. 4: (Color online) a) Left-continuous process Uj{t); b) Com- 
pound counting process Z{t) = Uj[S^{t)] as a new arrow of time; 
c) ID trajectory of standard Brownian motion (parent process); d) 
ID trajectory of Brownian motion subordinated by Zuit). Here ev- 
erywhere the index 7 equals to 0.9. 
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FIG. 5: (Color online) 2D trajectories of standard Brownian motion 
(a) and Brownian motion under compound subordination with 7 = 
0.7 (b). 



where the parent process X{t) is just the diffusion front cor- 
responding to the simple random walk R{t) given by (|7]), 
while the undershooting subordinator Z{t) coiTesponds to the 
limit of the rescaled counting process U^^ {v{t / Tf))) / g{Tf)) as 
''0^0 (with renormalizing function .g(To) chosen appropri- 
ately). Both the parent process and the undershooting sub- 
ordinator are well defined if the distributions of the spatial 
steps Ri and cluster sizes Mj satisfy some conditions, refer- 
ring to their asymptotic behaviors. In particular, if {Mj) < oo 
(or equivalently (fj) < oo) we have Z{t) = t/At. On the 
other hand, taking into account clustering with a heavy-tailed 
cluster-size distribution 

Prob(A/, > m) - (m/c)~^ 

m^oo 

with the tail exponent < 7 < 1 and some scaling constant 
c > 0, we obtain the compound form Z{t) ~ [S^{t/ At)], 
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where U~{t) = MvUx^t- U^{x) is the left limit of the 7- 
stable subordinator U^{t), and S^{t) = \ni{t > : U{t) > 
t} is its inverse process. For a fixed t > 0, the compound 
undershooting subordinator U~[Sj{t/ At)] has the general- 
ized arcsine distribution rescaled by t/ At. The corresponding 
probability density function reads 

p Jt, t) = T^-\t/At - t)-^ , < r < t/At . 

TT 

(10) 

It is easy to check that all moments of the random variable 
U-[S^{t/At)] are finite. 

The numerical approximation of the process [S^ [t / At)] 
is shown in Fig.|4] For simplicity, we take At ~ 1. To simu- 
late the process U~ [S^ {t)] , one only needs to generate the val- 
ues [/^(nAr) (see Fig.|4ja)), where n = 1, 2, . . ., and At is 
the step length. This can be carried out by the standard method 
of summing up the independent and stationary increments of 
the Levy process llzsll . Then the approximation of the pro- 
cess U:^[S~f{nAt)] is a simple continuous-time random walk 
in which each waiting time is exactly equal to the jump. In 
the result, Fig. H^b) demonstrates the process [S^ {t)] non- 
decreasing in time. 

As far as the parent process is concerned, for simplicity, 
but without loss of generality, we can consider the standard 
Brownian motion B{t) that can easily be obtained if we take 
into account the spatial jumps satisfying conditions = 
and (Rf) < 00 (see Fig.Hfc)). Then the subordinated process 
B{U~[S-f{t)]) represents the corresponding diffusion front. 
Its sample paths are shown in Fig. HJd). The main feature of 
the subordinated process consists in random jumping along 
random trajectories of the parent process. For better visual- 
ization of this feature we can consider a walk on a plane with 
2D vector jumps satisfying analogous conditions that leads to 
2D Brownian motion as a parent process, see Fig.|5] The dif- 
ference between Figs. [Sja) and (b) is that the walker, moving 
along a Brownian trajectory in presence of the compound sub- 
ordination, stops from time-to-time and overjumps through 
intermediate positions in the Brownian trajectory. In other 
words, the trajectory of Fig. |5lb) is something like a parti- 
tion of the trajectory in Fig. [Sja) on random intervals. The 
appearance of the overjumps, as applied to the relaxation of 
semiconductors with metastable defects described in Sec. |II] 
is connected with long-range interactions among DX centers. 

The probability density function of B{U^[S.y{t/ At)]) 
reads 



p{x,t) 



where p^{x, t) = 



p^{x,T)p^{t,T)dT , 



(11) 



Taking the inverse Fourier 



transform with respect to x for p(a;, t) in Eq.(fTTI). we derive 
the relaxation function = (^^-'^kR" {t}^ that characterizes 
the temporal decay of a given macroscopic mode fc. We get 



0(0 



sm 7r7 



t/At 



-(feV2)r^7-i(i/At-r)-^dT, 

(12) 



where fc^/2 may be denoted as a characteristic frequency ujp 
of the relaxing system. The integral representation of the re- 
laxation function may be substituted by a series expansion ob- 
serving that the relaxation function ( fT2] i can be expressed as 

(j>{t) = Ej ^{—ujpt), where 



fe=0 



{l.k)x^ 
T{ka + P)k\ 



a,/3 > 0, 



is the generalized Mittag-Leffler function |26ll. Here (7, fc) = 
7(7 + 1)(7 + 2) . . . (7 + fc — 1) is the Appell's symbol with 
(7,0) = 1,7 7^0. 

Based on the Laplace image of (t>{t) with respect to t we 
find the kinetic equation of relaxation in the pseudodifferential 
form 



dt 



r(i-7)' 



with the initial condition 0(0) = 1. For the experimental stud- 
ies of the dielectric spectroscopy data (dielectric permittiv- 
ity or the corresponding susceptibility) the frequency-domain 
representation of the latter function 



(/3*(w) = / e 
/o 



M}\ dt 

dt J 



(13) 



is of interest. Taking the Fourier transform, we get 



^*{lo) 



1 + ioj/ujp 



, 0<7< 1 



(14) 



Let us note that formula ( fT4b coincides with (O for a = 1. 
Moreover, for 7 = 1 it takes the form of the Debye re- 
sponse, corresponding to the spatio-temporal clustering with 
(Mj) < 00, where the undershooting subordinator becomes 
deterministic: Z{t) = t/At. 

Considering the dielectric susceptibility 
X* (uj) = x' {^) - ix" {^) = x(0)</3*(w) with (p*{u) as 
in (fl4] i one can show that the real and imaginary parts of 
the susceptibility fulfill the following frequency-independent 
relations: 



lim , , , , 

x'(o)-x'(w) 
x"i^) 



tan {^) , 



lim 



X'(^) 



(15) 



This means that for small lu the gallium-doped 
Cdo.Q3Mno.01Tc mixed crystals obey the energy crite- 
rion llSll . while for large uj the high-frequency energy lost per 
cycle does not have a constant relationship to the extra energy 
that can be stored by a static field. The similar feature takes 
place in the Cole-Davidson law where, however, the energy 
criterion is obeyed for large uj only. As it is well known lUSll . 
the presence of the universal law in the relaxation in gallium- 
doped Cdo.9gMno.01Te mixed crystals is caused strictly by 
many-component interactions. This effect is just accounted 
for using the compound counting random processes. 
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IV. CONCLUSIONS 

In the diffusive framework we have brought to Hght the 
fractional dynamics of gallium-doped Cdo.g9Mno.oiTe:Ga 
mixed crystals, for which the particular case of less typical 
relaxation behavior with the high-frequency power-law expo- 
nent 1 — n = 1 has been detected. We have proposed the sub- 
ordination scenario of anomalous diffusion underlying the ob- 
served pattern, based on random, heavy-tailed spatio-temporal 
clustering procedure applied to the random walk hidden be- 
hind the standard Brownian motion. Such a procedure has 
led to a substitution of the renewal counting process, used in 
the classical CTRW model, by the undershooting compound 
counting one. As a consequence, the explicit stochastic struc- 
ture of the diffusion front as a subordinated Brownian mo- 



tion has been obtained. The resulting frequency-domain re- 
laxation function has been shown to fit the studied dielectric 
spectroscopy data. 
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